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ABSTRACT 


Quasi -specular radar scatter from geologic surfaces displays a 
variable wavelength. A, dependence in apparent surface roughness, 

ranging between a ^ and o for .01 < A < 1 m. 

The strongest changes in a with wavelength are observed in lunar 
mare, while scatter from lunar highlands and most of Mars' equatorial 
region is wavelength independent. Commonly used, gently-undulating 
surface models for electromagnetic scatter predict no wavelength depen- 
dence. Wavelength dependence will occur whenever a significant fraction 
of the surface has local radii of curvature comparable to the observing 
wavelength. This condition can be determined by comparison of the value 
of the integrated surface curvature spectrum with the radar wavenumber, 
multiplied by a constant that depends on the geometry. 

Variations in curvature statistics calculated from photo gramme trie 
reduction of lunar images are consistent with the observed variations 
in quasi -specular scatter at the same locations. Variations in the 
strength of the wavelength dependence are correlated with the sizes 
of lunar craters that lie near the upper size limit for the local steady 
state distribution. This correlation is also consistent with variations 
in the curvature spectrum calculated from crater size-frequency distributions. 
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Introduction 

Earth-based radar observations of the moon and terrestrial planets 
are dominated by a small area of intense scattering in the center of the 
planetary disk ( Evans , 1969). Similar observations that use space probes 
to achieve oblique-scattering, or bistatic, geometries give equivalent results 
in which the strong scattering area is always approximately centered 
on the point where theanglesof incidence and reflection are equal with 
respect to the mean surface ( Tyler and Howard , 1973). Earth-based ob- 
servations represent the special case of normal incidence and reflection. 
The polarized scatter from , the dominant bright area is termed quasi- 

specular . 


The scatter from portions 

of the surface that do not approximate the specuTar conditions is observed 
to be insensitive to the detailed geometry, and 

is termed diffuse . Diffuse scatter generally obtains in radar mapping 
from the earth where backscatter at moderate to high angles of incidence 
is -employed. Quasi-specular scatter is dominant in studies of scat- 
tering laws for' the moon and planets, in oblique scatter 
experiments, and in radar studies of Mars, where fundamental sampling 
considerations limit earth-based analyses to a small area surrounding 
the sub-radar point. Only quasi-specular scatter is considered in this 
paper. - 

Quasi-specular scatter can be modeled by random surface models which 
have considerable roughness on scales that are much larger than the radar 
wavelength. That is, both the horizontal and vertical roughness scales 
are many wavelengths in extent. Sub-surface reflections, shadowing, and 
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edge effects can be shown to be unimportant for many practical cases, 
and are neglected. Typically, these surfaces are characterized as 
gently undulating. The principal contribution to the scattering from such 
surfaces arises from those points on the surface where the local slopes 
are properly oriented to produce a specular ray path — angle of in- 
cidence equals angle of reflection — between the radar transmitter and 
receiver, whether these are co-1ocated or separate. The assumption that 
the surface roughness scales are large with respect to the radar wavelength 
leads to the analytical result that scatter from gently undulating 
surfaces should be wavelength independent. 

Quasi -specular radiov/ave scattering from lunar and planetary surfaces is 
variable in its behavior with wavelength. Quasi-specular scatter from lunar 
plains units consistently reveals a strong wavelength dependence in apparent 
surface roughness over the range of centimeter to meter wavelengths, while 
most other lunar units show little or no variation in quasi-specular scatter 
with wavelengths in this range ( Tyler and Howard , 1973). Measurements of lunar 
samples show the electrical properties of the rocks to be nearly constant over 
this wavelength range. With a few exceptions, radiowave scatter from the 
equatorial region of Mars also appears to be independent of wavelength within 
this range ( Downs , et al., 1975; Downs et al . , 1978; Simpson et al . , 1977). 
Except for the observed variation with wavelength, the quasi-specular scatter 
from lunar plains is indistinguishable from other lunar units, or Mars, or from 
expectations based on theory. These other surfaces do exhibit somewhat differ- 
ent radar slope frequency distributions, however. The variations in radar 
roughness with wavelength correlate well with similar variations in photo- 
grammetric roughness with sampling scale. 
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Wavelength variability has been discussed generally in terms of a filtering 
property of the scattering process wherein the quasi-specular scatter is 
unaffected by surface structure which is small with respect to a wavelength 
( Hagfors , 1966). We have attempted in a previous paper ( Tyler , 1976, 
afterwards called Paper I) to give quantitative definition to this 
process by estimating the conditions under which a real surface is 
statistically equivalent to the type of surface assumed in the models 
of quasi-specular theory. Real surfaces that are well approximated in 
this sense are expected to display little or no wavelength dependence 
in radar properties, while those that violate the modeling assumptions 
are expected to display some variation in radar characteristics with 
wavelength. Very little other work in this area has been carried 
out. 


In the remainder of this paper v/e discuss a comparison between the 
observed viavelength dependence in lunar radar scatter data and statisti- 
cal measures of the surface derived from orbital images. Two types of 
comparison are possible. First, at a limited number of locations we 
are able, through the use of photogrammetrical ly determined heights, to 
calculate the distributions of height variance with inverse lateral 
surface scales, or frequencies. The surface roughness parameters that 
are important in determining wavelength dependence are calculated from 
these distributions. Results of these calculations are compared with 
oblique-scatter bi static-radar results from the same location on the 
moon, or from within a contiguous, independently mapped geologic unit. 
This comparison suffers primarily from inadequate surface resolution 
in the orbital images which leads- to an extrapolation of the 
height- variance spectra to small scales. The second method is based 
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upon a statistical model of lunar surface erosion that is used to determine 
relative surface ages. The frequency-size distributions of lunar craters 
in plains units are well behaved, and characterized by a particular crater 
diameter, C^, that separates the members of the "old", and "new" or 
"steady state," crater populations. This crater diameter is readily de- 
termined from images. The observed crater populations can be related to 
the height variance spectrum theoretically through the model. 

Again, the height-variance spectrum can be related to the expected 
behavior with wavelength. Bi static-radar data are used for the comparison. 
This method suffers from area sampling problems associated with the 
determination of C^, and the requirement for other assumptions that have 
not been tested. However, the model on which it is based has been 
validated to lateral scales of 10 cm or less by direct observations of 
the surface at Surveyor and Apollo landing sites, and thus is applicable 
directly to the surface scales of concern .here. 

More thorough discussions of quasi-specular scatter and its relationship 
to planetary radar data can be found in Paper I and 
elsewhere (e.g., Beckmann , 1963; Hagfors , 1966; Barrick , 1968, 1970; 

Evans , 1969). Radar procedures and results used here have been described 
by Tyler and Howard (1973), and by Tyler et &1 1973). Detailed comparisons 
of oblique-scatter radar data with other remote sensing data are also of 
interest in understanding the importance of various geologic surface 
expressions on radar scatter, and are given elsewhere ( Moore et al . 

1975, 1976). 
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Requirements for Radar Model 

Consider a two-dimensional surface, ?(x,y), generated by 
an isotropic gaussian random process with mean zero, = 0, and 
mean square height, =■ h^^ . Let represent the 

surface height correlation function evaluated at separation 
T = [(x-|-X 2 )^ + (y-|-y 2 )^]^'^^- The quantities 

( 0 ) = 


5 “ °^xx 

^xx^xx 


= 0 , 


yy 


are the variances U and , respectively. These quantities can 
be expressed as weighted integrals of the surface height-variance spectrum, 


0jf = 0jf = /° q^ S(q)dq 

*0 


= 0^/2 = (tan^y) /2 


*^xx f 

*0 


( 1 ) 


2 

where is the variance of the total slope or surface tilt, and 

the height-variance spectrum. 


S(q) = 


CO 

2tT J' R(t) J^(27TTq)TdT 
0 


and autocorrelation function 


( 2 ) 


R(t) 



0^{2irTq)qdT 
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3 5 

form a Hankel transform pair. VJs refer to q S(q) and q S(q) as 
slope and curvature spectra, respectively. 

Electromagnetic models for gently undulating surfaces always assume 
that the radius of curvature of the surface is everyvi^here much larger than 
the wavelength of interest. Paper I gives the probability distribution 
of surface curvature at the local specular points as. 


pci'll 2 I - '^12 ^ “ 

2.8[erfc(^) + .71 erfc(^y'^^ + .87 exp ("M) ]; ( 3 ) 

2 

where r-^^ 1 s the product of the principal radii of curvature, and 

Y is the tilt of the local surface at the specular point. This 

2 

function is plotted in Figure 1 . The probability that |r-j 2 j exceeds 

some particular value r.j 2 is a very rapidly decreasing function 

of p, for values of p less than about 1. The model will closely approximate the 

real surface, and vice versa, if P(|r^ 2 ^ [ < r^ 2 ^ ) is very small for 

an appropriate choice of threshold curvature, . 



7 


If the height variance spectrum is known to scales significantly 
less than a wavelength the value of P(-) can be estimated, 

O 

The quantity o is easily obtained from { 1 ), while cosy - 1 
for quasi-specular scatter. In our previous discussion we suggested 
that the Rayleigh criterion r-jg ^ k - Ig-I , where X is the radar 
wavelength, provides a firm lower bound on the radius of curvature. 

With these substitutions the parameter p becomes 


U 

% 



(4) 


which depends only on the wavelength and the curvature spectrum. 


2 

The infinite upper bound in the computation of a is not strictly re- 

AA 

quired, because it is known that roughness on very small scales with respect 

to a wavelength does not affect the specular scatter significantly ( Barrick , 

2 

1970). . Thus, o can be replaced with a less rigorously 

AA 

defined quantity 


provided that 


q^,<M 0/X 

a’^ = 12ir^ ( q^ S(q)dq, 

XX J 

0 



(5) 


where the integration is carried out over surface structure of scale 
greater than about X/10. Because parameter p always increases 
with X , every surface will violate the assumptions of a gently un- 

O 

dulating model for sufficiently large wavelength unless decreases 

XX 

_2 

at least as rapidly as X . 
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Quasi-specular scatter vn'll be wavelength independent over any 

2 

range of wavelengths for which y = — 9 about 

aA 


2tt" 


1 , for a fixed X . As an example, consider two wavelengths X-j 

I 

and X2 with X-j>X2 = X . Then if 


and 


y 


2,2 


y 


1,2 


2 ir' 


2 ^xx ^^^2^ 


27r" 


-XX ' (^2) 


( 6 ) 


are both less than about 1 (y^ 2 will be greater than y2 2^’ The surface 
can be well -approx-i mated -by Td gently undulating model and there will be no sig- 
nificant vavelength dependence^ over this range. Stated another 
way, given ■ a (X^) such that y„ „< 1 , the upper limit over X, 
over which little or no wavelength dependence is expected can be easily 

I 0 

found. Surfaces with sharply declining curvature spectra have a ^ (X) 

XX 

2 

that rapidly converge to their final value. ‘ - These surfaces display 

constant quasi-specular scatter for all wavelengths less than some 

. ' 2 

critical value. The rapid convergence of a corresponds to a lack 

XX 

of small scale surface structure. Conversely, fora real surface at any 

X' ^ '2 

wavelength such that y* = — ^ o (X‘ } > 1 , the assumptions of the 

2/ XX 

gently undulating model are not fulfilled and wavelength variations in 


the scatter can be expected. However, no general, rigorous analytical 
method is available to calculate these variations explicitly. 
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In Paper I we estimated the wavelength dependence of the 

apparent surface roughness, as sensed by radar, for several examples of 
power law height variance spectra under the assumption that a fixed 
value of u = u' defines an upper bound to the spatial frequencies that 
contributes to quasi -specular scatter. These results are summarized in 
Table I and will be referred to below. Note that there is a marked 
change in the strength of the wavelength dependence for values of the 
spectral index n between 3 and 4. For larger values of n, a is 

A 

virtually independent of wavelength. Also, there is only slight change 
in the power of the wavelength variation, from -1/2 to -1/3 as the 
height variance spectrum changes from a constant to an inverse cube law. 
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Comparison with Height- Variance Spectra 

The work described in this section was based on detailed lunar topography 
as determined by photo gramnietry, supplied by H. Moore and his colleagues at 
the U.S. Geological Survey. Their results are based upon Apollo photographs 
taken from lunar orbit. Detailed descriptions of the photogrammetric procedures 
and error sources have been given elsewhere ( Moore and Wu , 1973; Wu et al. , 1973). 

We obtained the data in numerical form as a sequence of planar location 
coordinates and associated heights. Samples were separated laterally by 
25 m, which was the closest spacing for which reliable estimates of slope 
could be obtained. Typically, the height was given for each of three deter- 
minations at each locating coordinate. The dispersion among these repeated 
samples was used to estimate the statistical sampling error and this error 
was converted to a value for the spectral noise floor under the assumption 
that the errors at different locations were uncorrelated. Several hundred to 
about one thousand samples were available at each sampling site. The sampling 
locations are arranged in systematic patterns that vary from site to site. A 
linear sampling traverse was most common, but rectangles, ’'L"s and zig-zags also 
occurred. These patterns were chosen to maintain all samples vn'thin a similar 
terrain type, or to avoid a feature that was judged by the geologists to be non- 
representative of the unit of interest, such as a bright crater ray in otherwise 
uniform mare material. Uniform area samples from which two-dimensional height- 
variance spectra can be computed directly were not available. 

Since our interest here is in the high frequency behavior of the 
height- variance spectrum, all data were treated as though they were obtained 
from a single linear traverse. In the few cases of folded patterns, this 
procedure distorts the low frequency portions of the height- variance 
estimate, but is not believed to introduce serious error for our purposes. The 
height-variance spectra computed from these data are detrended and have 
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30 degrees of freedom. Typical results are given in Figure 2. 

The asymptotic behavior of such spectra was estimated by two methods. 

Let S(q ) represent a one-dimensional height variance spectrum estimated 
X 

from the data. It was assumed ( on the basis of previous results) 
that the asymptotic behavior has the form 


S(qx) = 


(7) 


The constants C, t were' obtained first by least-mean -square fits to log (C) 

-t log (q ) between limits q ^ % > which were selected by 

X %i n ^max 

inspecti'on of plots such as those given in Figure 2. The fitting region 

was chosen to avoid low freouency turn-overs and, occasionally, erratic 

fluctuations in the extreme high frequency portions of the spectrum. Typically, 

q = .001 cycle/m, while q was chosen as the highest available frequency 
%i n ^max 

(.02 cycle/m) unless the estimates of sampling noise and the appear- 
ance of the spectrum both indicated a problem at the higher frequencies. 

As a check on this procedure we also determined the best linear fit 
to 


log (S{q^) - N) (8) 

where N represents the spectrally uniform contribution that would result 
from uncorrelated photograminetric- reading errors. The value of N for’ . 
which the fitting 
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error was minimuin was taken as an estimate of the noise floor. Corres- 

ponding values of C, t were taken to represent the asymptotic behavior of 

S(q ). In this case q . = .002 cycles/m and q^ = .02 cycles/m. 

^ ^ ^min Xmax 

The values N and the corresponding estimate obtained from the sampling error 
generally agree within a factor of about 2. 

As a final estimate of C, t the values of the spectral index t 
obtained by the two methods were averaged. The spread in the two values was 

used as an estimate of error. Because the constant Cwas believed to be much 
better determined than the slope, and the final results are not sensitive to 
its precise value, only the second determination of this quantity 
was used. 

To ,be useful in a two-dimensional power-law formulation of the 
surface roughness, these results must be interpreted in terms of the 
two-dimensional spectrum 


S(q) = k q“P 


( 9 ) 


It is shown in the Appendix that the relationship between the power law spectrum 
SCq^,) based on one-dimensional sampling and an underlying two-dimensional 
power-law S(q) is 

p = t + 1 

= .1087 


k 


2<t < 5 


00 ) 
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The relations (10) were used to express the^.spectral estimates S{q ) ob- 
tained from the linear fitting procedures discussed above in terms„.of the 
twp-dimensTonal parameters k,p. 

Table II gives the results for all but one of the sites for which 
oblique-scatter data are also available. The rejected site gives 
anomalous results for k,p. Examination of the orbital images shows 
that the photogrammetric results are contaminated by a fresh crater. 

Values of the error estimate for p are given in parentheses adjacent 

to the quantity. The indicators and a and'^S refer to one case for which 

only the second estimation procedure was used, and to several cases for 

which the two estimation procedures yield the same result. Note that the 

final values for the spectral index generally distribute themselves between 3 and 

about 4, and cluster near the extremes of this range. 

2 

Table II also gives the surface mean-square height, h^ , calculated 
directly from the original height sample data. Since k, p represent 
only asymptotic behavior, a comparison of the spectrum kq ^ and 
h^ is a measure of the low frequency surface behavior. 

For reasons that are discussed in the next section of this paper, 
it is not expected that the high frequency tails of S(q) will turn 
upwards, i.e., that p will decrease for frequencies higher than those 
sampled. If the value of p changes, it is expected to increase, i.e., 
the spectrum will remain the same or become steeper. Still, to our 
knowledge there are no data on the height-variance spectrum of the 

lunar surface at frequencies higher than those shown here, so extrapolation 
to much smaller scales than those sampled directly is required. 
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The parameterized asymptotic height spectra can now be 

used to determine the surface parameters of interest to radar scatter. 

2 

In particular, we could compute o ' (A) directly>and from it the para- 

meter p', to determine the consistency of the surfaces sampled with the assumptions 

of the quasi'specular model. However, for purposes of comparison with radar data it is 

more convenient to use S(q) to determine the spatial frequency q^ (A) 

2 

at which u' exceeds some threshold. Redefine o' in terms of 

XX 

specific values for the upper and lower bounds. 


'2 5 

Oxx = 127T 


q S(q) dq 


(11) 


v/here q^ is the effective upper bound and q^ is the effective low 
frequency cut-off required for a power-law spectrum. Substituting 
kq~P for S(q), 



( 12 ) 


For our data, q^ is computed by extrapolating kq“P to low frequencies 
until the integrated value of the height spectrum is equal to the mean- 
square height, i.e., by solving 
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^0 = J q’S(q) dq 

% 

for q^. It is easy to show that the result is insensitive to the 
upper limit for the spectra S(q), of interest here. Values of 

_c ‘ _5 

are typically 10" cycles/m, with a maximum value obtained of 3.5 x 10 
cycles/m for Mare Serenitatis in Table II, Values obtained for q.j are much 
larger. Thus, the second term in the expression for < can be neglected, 
with the result that 


qf = 


(6-p)o, 


,2 


XX 

■T2tt^-K 


6-p 


p < 6 


(13) 


The asterisk is used to denote the value of q.j estimated from the data, 
i.e., for the specific values of k,p obtained from reduced surface height 
sample data. Referring to Figure 1, we take p' = 1 as a reasonable 
upper bound for the conditions under which a real surface can be approxi- 
mated by a gently undulating model approximation. For this condition about 90 percent 
of the specular points fulfill the required approximations; for smaller 
values of y‘ the degree of approximation rapidly becomes much better. 

Further, we take cos% = .8 as a typical value for this quantity to 
account for components of the scatter somewhat removed from the centroid 
of the scatter. 
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Then, since 



? ' 

M-.- 

4 ,2 
cos Y ^ 


v;e have 

a^x^(13) = 910 (14) 

(116) = 12 

Expression { 13 ) above has been evaluated for the locations of interest 

' 2 

using these values of a at A = 13 and 116 cm, the two wavelengths 

XX 

of the Apollo bi static radar observations; the results for are given in 
Table II.. The range of corresponding to the spread in values of 

p is about a factor of 2; extreme values are given in parentheses for 
q-j*^ (13) only. A similar range of variation will result for (116). 
Typically, the values of q-j* lie between about 1 and 100. A notable 
exception occurs at the Apollo 17 landing site, where q| (13) and q| (115) 
are 550,000 and 11,000, respectively. This particular surface is extremely 
smooth at small scales. 

A second test for consistency between the gently undulating surface 
and the asymptotic height variance spectrum lies in the value of the calculated 
mean-square slope a for a particular upper frequency bound, q, , deter- 
mined by the curvature. Again, it is more convenient to compare the bounds 
for a particular case. It is easy to shov; that 
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p < 4 (15) 

„ 2 

C = iz£. _ X . . 

4tt (4-p) 

where the hat denotes the simultaneous solution with respect to curvature, 
a , and slope, o . That is, for a particular spectral index p, q, 

XX X I 

and k-j are the bound and constant that result in particular values of 
and * The quantities q-j , have been evaluated and 

are given also in Table II for particular values a = 12, corresponding 

XX 

2 

to y‘ = 1, A = 116 cm, and o = 0.02, v/hich represents a typical lunar 

X 

highland surface tilt of 11® rms. 

Values of q-j , k-j for other choices of these parameters can be easily 
obtained from those tabulated in Table II and the formulas above. 

The criteria for wavelength independence developed above under 
“Requirements for Radar Model" can be applied to the q-j* (116). According 
to our previous discussion, no wavelength dependence is expected for 
quasi-specular scatter if y-j 2 < 1 » or what is equivalent, if 
q* (116) > 8 cycles/m. That is, that the surface characterized to scales 
of 13 cm must also be gently undulating on the scale of 116 cm. This cri- 
terion is clearly satisfied by the surfaces at four locations in Table II 
and violated at seven others. There is one borderline case, Descartes. 
Generally, these surfaces are separated into two categories for which 
q* (116) 'V 2 and q* (116) 'v 10. 

/*v 

This separation is stronger when values of k, k-j and q-|*(116), 
q* are compared. Those surfaces for which y' < 1 are also those for 
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which the k-j , q-j required to bring the slopes and curvature simul- 
taneously into agreement are near the k, q^ (116) derived from the surface 
height data. For these sites k ^ k-j , although the k-^ are 
systematically too small by a factor of about 2, while q|" (116) 'x, q-j , 

/S. 

but the q-j may be systematically too large. These small systematic 

disagreements are not surpri si ag. Thefe are large differences between 

these pairs of values for the other locations. This is particularly 

evident in the differences between k and k-j , where k typically 
✓\ 

exceeds k-j by an order of magnitude or more. Most of the sites 
in this second category are in areas of low to moderate slopes. The 
discrepancies would be considerably larger if smaller values of a 
corresponding to the observed slopes for these locations had been used 
to estimate , k rather than the value = .02 actually used. Descartes 
apparently should be included with the first group of sites. 

Finally, Table II also includes values of rms slope obtained from scatter 
observations with Apollos 14 and 15. In simplest terms, the values 0^(0 given 
are direct measures of the angular one-half power width of the obliquely scat- 
tered echo, corrected for certain geometrical factors, and reduced to a unidi- 
rectional rms slope on xhe basis of a gently undulating surface model. The para- 
meter, 13 or 115, indicates the wavelength of observation in cm. Complete descriptions 
of this experiment, the underlying theory, and the data reduction process have 
been given elsev/here ( Tyler and Howard , 1973; Tyl er et al . , 1973). 

Some further comments are needed for several of the sites. There are no 
radar data of this type for the Apollo 17 landing site; this site 
v/as included in the table only as an example of an extremely smooth 
surface sampled by photogrammetry. There are two sites for which 



19 

no radar data were directly available, but where geologic maps 
(Eggleton, 1965; Howard and Masursky , 1968} 

indicate that the radar has sampled the same unit within a short 

distance of the photo grammetric traverse., These sites are 

Cratered Terrain and Mare Cognitum. Mare Cognitum is a large, uniform 

area with little complexity. Cratered Terrain is embedded in highland 

materials. One location, Cayley, is included although radar data are 

given only for A = 13 cm. This unit is included because it was -observed 

directly by Apollo 16. However, there were severe experimental problems 

with the 116 cm wavelength data from the Apollo 16 experiment, so that no reliabl 

results are available for a (116). The Apollo 16 data show no significant 

X 

change in the scatter at either 13 or 116 cm wavelength between the 
Cayley unit which is located in the crater Ptolemaeus and the highlands 
material surrounding it. Since highlands units observed with both 
Apollos 14 and 15 generally display wavelength independent scatter, 
this Apollo 15 observation is taken as evidence that the scatter -from Ca^^ley is also wave- 
length independent. Further, the values of o (13) for Apollo 16 are 
known to be systematically low with respect to those of Apollos 14 and 
15, but the cause of this problem is not known. 

All locations for which the photogrametri cal ly sampled surface is 
apparently consistent with the gently undulating surface display v/avelength 
independent quasi-specular scatter, as expected. Of the other locations, 
only Cratered Terrain and Mare Undarium fail to show a strong wavelength 
dependence, although a weak variation is present. The remaining sites 
typically have a 1.5 to 2.0 to 1 variation in a between 13 and 116 cm 

X 

wavelength. On the basis of our previous work, a 2 to 1 variation 
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would be expected over this range for an inverse cube law dependence 
in S(q) (v. Table I). 

These results show a general quantitative agreement between the 
onset of wavelength dependence in radar observations and the consistency 
of the surface with the radar model as determined from 
independently measured surface height statisics. If 
the assumption that the variance spectrum of the lunar 

surface only becomes steeper at smaller scales is correct, then the 
unexpectedly weak variations in a at Mare Undarum and Cratered Terrain 

/V 

locations can be explained also. . : 

Unfortunately, Table II exh'austs the available sets of co-located 
lunar topographic data and radar observations. There are no appro- 
priate lunar surface data of this type other than those considered here. 
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Comparison vnth Values 

One of the more striking characteristics of the lunar surface 
is the ubiquity of craters and the absence in orbital photographs 
below a scale of about 1 km of any apparent change in the distribution 
of the crater population. This latter characteristic was predicted by Moore 
(1964) and others ( Shoemaker , 1965). It is generally believed to result 
from a saturation bombardment of the surface by meteorites, to the level 
where old carters- are destroyed as fast as new ones are created. The 
size of the largest crater that belongs to this steady population in- 
creases with time. Craters larger than this limiting size have a 
distinctly different distribution whose origin is not well understood. 

The size- frequency distribution of lunar craters can be determined 
from photographs. Figure 3 gives a schematic 

illustration of the cummulative crater population for a typical surface. 

-2 3 

The region represents the steady state condition, while the C 

is typical for craters greater than a particular size, called C^. For 

sizes greater than about 1 km, there is a second break to a third 

slope which is not of interest here. 

Lunar crater size-frequency distributions have been studied exten- 
sively, and are well represented by the curves just described (v., e.g., 
Soderblom and Lebofsky , 1972; Gault , 1970), As indi- 

cated C^ increases with time, and the breaking point moves along the 
extrapolated line. The underlying theory for interpretation of 

these distributions has been confirmed by numerical modeling and experi- 
mentally in the laboratory so that the functional form is well established. 
Of considerable importance here, the 5 portion of the distribution 
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has been verified to scales at least as small as 10 cm by in-situ observations 
of the lunar surface (Morris and Shoemaker » 1970). This observation is 
thought to be representative of plains surfaces on the moon, but may break 
down in extremely rough areas where other mechanisms, such as downslope trans- 
port, also modify the surface at small scales ( H. J. lioore. private communi- 
cation, 1977).' Because a determination of requires careful counting 
of craters, a second measure of surface age based on recognition of the 
eroded morphology of large craters has been developed ( Soderblom and Lebofsky, 
1977). This measure is Dj^, which is defined as the diameter of the largest 
crater that has eroded to an internal slope of 1°. The quantity Dj^ is 
readily determined from lunar photographs by identification and measurement 
of certain craters with a particular shadow geometry, with a simple correction 
for sun angle. Values of are available for most of the Apollo 14 and 
15 oblique-scatter data from plains units ( Moore et al .. 1975, 1976)^ The 
quantities and are in approximately constant ratio, = 1.7 
( H. J. Moore, private communication, 1977). Because of the simple relationship 
between and the size-frequency distribution, it is easy to derive an 
approximate relationship between Dj^ -and the form of the height-variance spectrum. 

We assume that all craters of the same size give rise to the same spectral 
contributions to the overall height- variance spectrum, i.e., that on the aver- 
age the height-variance spectra of individual craters of the same size are the 
same, and that the spectra of individual craters scale in accordance with the 
variation of a crater of variable size but fixed shape. Within the steady-state 
fraction of the distribution this is clearly not the case, because craters of 
any particular size are found in all stages of erosion. However, this factor 
will tend to smooth the small scale surface more rapidly than would be 


the case under our assumption. 
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Let f(p) be an assumed crater profile, where p is the radial 
coordinate, then the crater scales as af(p/a) if both vertical and hori- 
zontal dimensions increase by the same factor a. The Hankel transform of 
f(p) is 


H[f(p)] = s(q) = 27 t J p f(p) JQ(27rpq) dp (16) 

0 

from which it is easy to show that 


H[a f{p/a)] = s(aq) 


07 ) 


_ •p 

The two-dimensional height variance spectrum, S(q) = js{q)i , then 
goes as a S(aq), where the bar denotes the spectrum of an Individual 
crater. Therefore, for a crater of size the Individual contribution 

r _ 

to the spectrum scales as ? S(^q) relative to the spectrum of an in- 
dividual crater of size unity. 

If all craters are statistically independent, then the spectrum for 
the total surface is 


S(q) = J P(0 i S(5q) (18) 

0 


where p(5) is the marginal probability density function per unit area that 
a crater of size C,5+d? is present. This simple form results from the 
assumption of statistical independence, so that the 
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spectra of the individual craters add incoherently. Other assumptions 
regarding the relationships among the various individual craters or other 
scaling laws lead to drastically more complex forms. This topic deserves 
further investigation. 

Describe the cumulative size-frequency distribution as 


P(C) = a 


5 < Dr 


a D ' 


Dq < ? < D^ = 0$ = Dj_/1 . 7 


(19) 


a D 2 D ■ 


D|‘ < C 


Where D^‘ is introduced as a convenience. Then, the marginal density 
function is 

P(C) = :^0“‘P(C)) = 0 ? <'Dq 

t 

2a D ^ 

= (20) 

3a D 2 D^' 


D|^ < £ 
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For this model a-1 > so must be interpreted as the size of the 
smallest crater. 

We also assume a power law for S(q) 


S(q) = b. 


q < qb 


= b(qj^/q)^ ,q|^ ^ q 


(21) 


S(5q) - b, ^q < 


= b 



% ^ 


The frequency q^/^ here plays the role of a cutoff for the spectrum 
of the individual crater. Substituting in (17) to find the surface spectrum 


Hin(D^',q|^/q) 




S(q) = Cyj bC^ d? + C^- f b ?^"^(qj^/q)^ d^ 


D. 


■/ 

Min(Dj_,qb/q) 


Max(D',qb/q) 

+ C2 *jf b dC + ^2‘J b (qb/q)'" d? 

\ Max(D|_' ,qb/q) 


( 22 ) 
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? 2 

Where C-j = 2 , and = 3 D|^'. Clearly only three of these 

terms can be non-zero at any one time; which three depends 

on the relative sizes of D^‘ and qj^/q. Evaluating the integrals. 




3-v [5-x» 


(lim - (D^' j. 


< q 


V 5>^ 3,4 


(23) 


= T { '•’l T - ('’l 


3-V 




V f 3,4 



27 


For V > 4, S(q) q"^, ■ ^ 

For large values of v, for which S(q) is approximated by 

{ 

I 

S(q) = b, ^ ^ ^b 

= 0, qj^ < q 

a more simple form is obtained* ^ 



Clearly, for all cases for which S(q) unconditionally converges 

i.e., V > 4, S(q) is of the form of a two-segment power law with 
indices 3,4 over most of the range of spatial frequency q. This char- 
acteristic arises from the form of the crater size-frequency distribution 

"2 "3 “ - 

(i.e., ? , ? ) and the edge, or break in S(q). The simpliPied result 

corresponding to large values of v is plotted in Figure 4. 

The critical parameter in this result is the ratio of the break 
frequency qj^ to Dj^’. There are no data on the variance spectra of 
individual lunar craters ( H.J. Moore , private communication, 1977). But, 
a small nuclear test crater, "Danny Boy", has been sampled adequately to demon- 
strate significant roughness to spatial frequencies that extend to many times 
the inverse crater diameter (Moore et al . , 1974). This characteristic 


is also readily apparent in images. Examination 
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of lunar images suggests that fresh craters have significant roughness 
to perhaps 100 times the inverse crater diameter (v. e.g. Moore , 1971). 

It is clear from experiments that the details of a particular cratering 
event, especially the nature of the target material, are of paramount 
importance in this regard (v. e.g. Oberherk . 1975). 

For our normalized spectrum S(q), the value of the break point is 
just the highest signficant frequency measured in inverse diameters. 

Typical values for Dj^' in lunar plains units range between about 100 
and 500. If the model for S{q) proposed here is approximately correct, 
then the break point for the spectrum will occur for q 'v 1 . This value 
is of the same order as the inverse of the radar wavelengths , , for the 

scatter data considered here. 

Note that the high frequency, portion of S(q) is con- 

trolled by the saturated, or steady state, region of the crater size- 

— - 

frequency distribution, while the low frequency segment, q , arises 

from the craters of size greater than C^. This occurs because the 

lav/ observed in the steady state represents relatively fewer small craters 

with respect to the total population than does the segment of the 

size-frequency distribution. The simplified spectrum for the cutoff 

S(q) goes to zero at q = which represents the highest frequency 

present in the crater population. Just below this maximum frequency, 

S(q) falls even more rapidly than q The mixture of crater ages in 

the steady state distribution will reduce the roughness of an average 

crater of a particular size relative to the scaled roughness of a crater 
-3 

in the ? population, so that qj^ should not be constant for the 
two portions of the population. The principal effect of this smoothing 
would be to displace the knee of S{q) towards smaller q. 
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As shown in Table I, we have previously estimated that 
-*3 "1/3 

S(q) q gives nse to X ' variation in a fX), while 

A 

S(q) -V q"^ leads to [to(x)]^^^> or a very weak wavelength dependence. 

Hence, radar scatter from surfaces of varying should display a 
variation in o (X) between these two 1 i mi ts, or possibly betvyeen 

A 

I/O 

X ' and X , as the changing places 

the observing wavelengths on different portions of the S{q) curve. 

Radar scatter and are compared in Figure 5 , where results 
for Apollos 14 and 15 are given separately. For each plot, the ratio of 
the rms surface roughness observed at the two Apollo bi static- radar 
wavelengths, p = o (116)/a (13), is plotted vs D. obtained from 

XX L 

Apollo and Lunar Orbiter photography. The values used are a slightly revised 
set of those previously published by Moore et al . (1975, 1976), based on more 
recent work. The data 'are restricted to lunar plains units. Each point represents an 
average value of D. and a over about 1° of surface arc along the bistatic 

L* A 

radar track on the lunar surface. A large area sample is required for 
reliable estimates of D^. In most cases there is very little or no 
variation in c within these areas. There is a clear systematic 
difference in the Apollo Hand 15 results( Moore et al .,1976). This difference 
has been corrected here by multiplicationof theApollo 15 o ratiobyl.22. 

A 

The Apollo 14 data show a clear, strong trend in p for 
between about 150 and 400 m; the Apollo 15 data show a similar result 
over the same range of D, values, but the relationship between p and o, 

^ L 

is much more scattered even though the formal errors in the data are about the same. 
Both data sets show p .5 for 150 m, and both data sets contain 

points for Dj^ > 400 m v/here p 'v 1 . There is a clear variation of p with’ 
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most of which occurs over the range 100 < < 400m. Further, 

the variation is of the correct sense with D^, and the minimum observed 

value of p'^».5 is approximately the correct value (13/116)^'^^= .48) , 

to be explained by a variation in S(q) of the type modeled here. In the 

Apollo 14 data, the minimum value of the ratio is p = .40, for Apollo 15 the 

minimum value is p = .46 (or .46/1.22 = .38 before scale factor 

“1/3 

correction). This seems to be in good agreement with X when it 
is remembered that these values of p are noisy. 

Linear regression v/as used to obtain an objective measure of the 
variation of p with D^. Other functional forms were also tried, 
but the data were not thought to warrant much more elaborate treatment. 

A value of 400m was selected as an upper bound to the fitting region 

on the basis of the Apollo 14 data , which have an apparent knee at 

about that point. The results are given in Table III for several 
combinations of data. 
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The two cases labelled Apollo 14 and Apollo 15 are the individual results 

for the complete data sets from those experiments for D|^ < 400 m. The 

increased noisiness of the Apollo 15 results appears in the much lower 

2 

regression coefficient, r , calculated for that case. Deleting the three 
points farthest from the regression line, the case labelled Apollo 15, 
except 3 points, increases the regression coefficient by a factor of 1.6. 

The fourth case, labelled Apollo 14 and Apollo 15, includes all data 
points and represents a reasonable fit of the regression line. 

The origin of the systematic differences between the Apollo 14 
and Apollo 15 data is not known. These differences are apparent in the 
radar data alone {cf. Apollo 14 and Apollo 15 results in Moore et a1 , 

1975, 1976). They may be related to a change in the experimental 
conditions, specifically the spacecraft antenna configuration, between 
the two sets of observations, or they may reflect real systematic differences 
between the Maria Serentitates and Imbrium, the primary source of Apollo 
15 data, and the Oceanus Procellarum, the primary source of the Apollo 
14 data. A third possibility lies in the sources of the images that were 
used to determine D^. The Apollo 14 values are based almost 
entirely on Lunar Orbiter results, while Apollo 15 values are 
based on those from Apollo ( Moore et al , 1975, 1976). Regardless 
of the differences, both sets of data display the same behavior, and the 
results from the combined set, Apollo 14 and Apollo 15, are taken as the 
best estimate of the variation, even though the regression coefficient 
is higher than for the Apollo 14 case alone. The changing wavelength 
dependence of radar scatter with lunar surface characteristics has also 
been presented elsewhere in forms that are readily compared with images 
and maps ( Tyler and Howard , 1973; Howard and Tyler , 1972; Moore et al . , 1975, 
1976). 
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Conclusions and Remarks 

The wavelength dependence of quasi -specular scatter is readily measured 
but has not been well understood. Scatter from most of equatorial Mars shows 
essentially no variation in the apparent surface roughness over wavelengths 
from a few centimeters to meters. Similar scatter from the moon shows little or no 
variation over the same wavelength range in highland areas, but does have 
an approximately 2:1 variation in most mare areas. It is believed 

that these differences in scatter are related to differences in surface 
roughness on approximately wavelength- si zed scales, but almost no quantita- 
tive work has been carried out. Paper I 

suggested a quantitative procedure to estimate wavelength dependence 

in terms of small scale surface roughness. as expressed 

through the height variance spectrum. 

0 ) 

This work has been an attempt/ to test the consistency of the gently 
undulating surface model for radar scatter and the natural lunar surface — 
which has not been done before — and to compare the results with observed 
scattering behavior, and (2) to verify the previous estimates of wavelength 
dependence given in Paper I. 

The circumstantial evidence for the dependence of the apparent 
roughness a (X) on the height variance spectrum S(q) is strong. 

A 

Expected radar characteristics based on estimates of the height variance 

spectrum are in good agreement with observations of radar scatter 

at the same locations. These comparisons suffer from the need to 

extrapolate to small roughness scales; but, it was also shown that 
-4 

S(q)aq is the expected form for the high frequency portions of 
the variance spectrum. On this basis, extrapolation of S(q) with 
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index p near 4 v/ould be accurate. For smaller values of p, the 
spectrum might soften, but is not expected to turn upwards. The 
results are consistent with this model. There are examples 

of strong wavelength dependence, and others where there is little or no 
variation with wavelength. Values of the diameter of a crater eroded to 
1 “ internal slope, are clearly related to variations in radar roughness with 

wavelength. The comparison of radar scatter with’ values is completely 
dependent on the'model for surface aging as a link' to S(q). While this 
model is v/ell established, the calculation of S(q) from that model is not 
rigorous. But a reasonable choice for the free parameter also 

leads to approximate quantitative agreement between the observed radar 
behavior and its variation with D^. Further work is needed on this 
problem. The realization that D|^ can be directly related to surface 
roughness on much smaller scales is apparently new, ^ 

Two different approaches to the data, one based on detailed 
samples at a single location, but limited in surface resolution to about 
25 m, and a second based on the statistical morphology of lunar craters 
over much larger areas, but which uses a surface aging model valid to 
scales at least as small as 10 cm, give similar results. Prior 
to this, there has been no demonstration that either the magnitude 

or approximate asymptotic shape of S(q) would yield values of 

2 2 

and y that were compatable with the gently undulating 

X XX 

model, or the filtered surface model for wavelength dependence. Neither 
has there been any model which predicted the variation in v/avelength 
dependence with Dj^, based on the changing position of the knee in 
S(q). Better understanding of this problem will require a new body of 
data from which radar scattering can be compared with surface statistics 
based on detailed measurements at sub-wavelength scales. 
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Appendix 

Consider a statistically isotropic two-dimensional surface described 
by a power spectrum. 


S(q) 


oo 



0 


R(t) r (2frTq) dx 


where R(t) is the auto-correlation function and S(q) is assumed 
to be asymptotically 


S(q) = K q"P, q.| < q 


i.e., to follow a pov'/er law for q sufficiently large. Then the one- 
dimensional spectrum associated with a profile of the two-dimensional 
along any fixed direction is 


S(q,) 
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2 2 2 

where q = q + q,, » and we have aligned the sampling direction with 
A y 

q . Evaluating this expression 


S(q,) = 


•q 


2 k 
P-1 


CO 

/ 




(Hqy2 )P/2 


so that if S(q) is power law, S(q ) is also with spectral index 

A 

reduced by 1. The integral 


f 


— 76 

is approximated by 2.30 p~* to less than 2% error, over a small 
range in p, 3 < p < 6. 

These results lead directly to the forms given in expressions (10 ) 


of the text. 
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TABLE I. ASYMPTOTIC WAVELENGTH DEPENDENCE 
FOR POWER LAW SPECTRA 


n 

0^ ^ 


1 

0 

X 2 

‘ 

1 

3 

n. X“ 3 


i 

4 

% [const. - £n{X)] ^ 


a, 0 £ q ^ 

S(q) = 

aq " q“", q-, < q 
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TABLE II. COMPARISON OF PHOTOGRAMMETRICALLY DERIVED 
SURFACE PARAMETERS ANO RADAR ROUGHNESS 
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LONG 

LAT. 

NAME 


k 


p 

q,*(i3) 

qi*(ii6) 

^1 

CM 

0 

q,(i3) 


% 

PHOTO 

10 

S/C SOURCE 
FOR 0, 


51®15'E 

Il^AG'N 

Mare Undarum 

100,000 

4.0x10"^ 


3.0(+.l) 

(TK 

5.714.6) 

1.4 

10 (9 8) 

1.6x10"®^ 

[]:$) 

5.5' 

4.5' 

2 5x10"^ 

9^89 ■ 

9484 

A-15 

5 

30°46'£ 

20“10’N 

A- 17 Landing Site 

110 

1.3x10"^ 


4.9(+ 1) 

5.5x10® 

1,1 xio'’ 

Y 

- 


- 

- 

2.6x10"® 

2755 

2750 

A-15 

21 

30°40'E 

20°15'N 

'Nest of A-17 

70 

3.4x10'^ 


3.9(B) 

99 

13 

27 

1,2x10'® 


5' 

5' 

1.6x10"® 

2755 

2750 

A-15 

25 

28“05'E 

20’00'N 

Mare Serenitatis 
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6.1x10'^ 

(Vb^) 

3 0(±.3) 

5 o(2'.'8) 

1.2 

10 (oiy 

1.6x10"®^ 

!Z li 
11.2) 

40 

2' 

3.5x10"® 

9565 

9560 

A-15 

3 

27“30'E 

20°10'N 

Mare Serenitatis 

noo 

9.2x10"'’' 
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11 Q!.) 
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4' 
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24025 'E 

10°40'S 

Northwest Flank 
of Theopinlus 
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4.3X10"'’ 


3.5(a) 

18 

3.2 

13 

2.2x10"® 


8' 

5' 

2.8x10"® 

10506/„(-, 

10504'^®' 

A-14 

15 

22“30'E 

tO'’20‘S 

Kant Plateau 
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4.3x10"^' 

ftf) 

3.0 (f.2) 

/8.81 
5 6(3.7) 

1.3 

10 (9 .^5) 

1.6x10'® 


7° 

' 5° 
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9“30‘S 

South of 
Zollner 

340 

3.1x10"® 


3,7(B] 

71 

11 

16 

2.1x10"® 


7° 

7° 

3.4x10"^ 

10624, 

10622(*®' 
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17 
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8°55'S 

Descartes 
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5.9x10"®' 

fid’) 

3 6(1.2) 

nool 
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7,7 

/19\ 
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7° 

7' 
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4563 

4558 

A-14 

20 
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Cratered Terrain 
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10 
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10 
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A-14 
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14 

19 
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4657 

4652 
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8 
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_ 
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a Only a single spectral index estimate available. 

6 Two spectral index estimates agree to two significant figures 

Y This result Is given as an example of a surface that is extremely smooth on a small scale. Values are not 
compared with radar because of limited region to which they apply, 

d Radar data from Apollo 14 at 2®30'E, 6°S, similar terrain. 

e Radar data from Apollo 16, 116 cm not used -values known to be in error, 

Cayley material in Ptolemaeus, Apollo 16 13 cm shows no significant change from heavily cratered material 
outside crater {USGS lA #80). 

1 Radar data from Apollo 14 at 24°W, 8°S from same map unit (IPH) as photogrammetric data, 
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TABLE III. 


p = a + bD|^, 


Case 

a 

Apollo 14 

0.06 

Apollo 15 

0.32 

Apollo 15, ex- 

cept 3 points 

0.29 

Apollo 14 and 

Apollo 15 

0.17 
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ANALYSIS OF p(Dj_) 
TOO < Dj^ < 400m 



0.0027 

0.87 

14 

0.0013 

0.28 

23 

0.0015 

0.44 

20 

0.0021 

0.57 
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Figure Captions 
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1. Probability that the magnitude of the product of the principal 

2 -^2 

radii of curvature jr^ 2 l less than parameter 

2 2 4 

r^ 2 ^os Y (see text). 

2 . Sample one-dimensional height variance spectra from lunar photo- 
gramme trie data. A. Data from Kant Plateau. B. Data from Cayley 
{see text. Table II for coordinates, parameters). Light lines 

are fits to curves according to expression (8) in text. The noise 
level for both curves is estimated to be 17 m /cycle/m ±10, under the 
assumption that the errors in the photogrammetric samples are uncor- 
related. The flattening of curve B for q >10 is believed due to 
noise from sampling errors. Neither curve is reliable for q 

X 

"3 

significantly less than 10 , due to the detrending procedures. 

3. Crater size-frequency distribution model for the lunar surface. 

Two segments of the distribution follow power laws with-indices 
-2, -3, and break point crater diameter C^. (This figure 
adopted from Moore, H. J., J. M. Boyce, and D. E. Meyer, "Lunar 
Impact Cratering in Recent Times, unpublished, 1977.) 

4. Height variance spectrum of model cratered surface. S(n) 
calculated from expression (23) in the text. Discontinuity in 
slope at n == 1 is real. Light lines give asymptotes, slopes 
are -3, -4. 

5. Comparison of surface rms slopes inferred from radar at 13 and 116 

cm wavelength with D^^, the diameter of the largest crater eroded to a 
1® internal slope. Data are from the lunar plains units only. 

Light lines give the results of fits to data for Dj^ < 400 m 
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Figure Captions (cont.) 

(see text, cf. Table III). Errors for are numerical average 
of formal errors for individual data points; errors in cr(116)/ 
Oy(13y are informal estimates. 

X 
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